Radiative transitions in charmonium from Nj = 2 twisted mass lattice QCD 
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I. INTRODUCTION 

Charmonium physics plays an important and unique 
role in our knowledge of Quantum Chromodynamics 
(QCD), which is believed to be the fundamental theory 
for strong interactions. In some sense, it is comparable 
to the hydrogen atom for atomic physics, the basic the- 
ory of which being Quantum Electrodynamics (QED). 
However, charmonium physics is much more involved in 
the sense that, due to its intermediate energy scale and 
the special features of QCD, both perturbative and non- 
perturbative physics are present. It is therefore an ideal 
testing ground for our understanding of QCD from both 
perturbative and non-pcrturbativc sides. 

Radiative transitions among various charmonium 
states are particularly important in the study of char- 
monium physics. Most charmonium ground states lie be- 
low the open-charm {DD) threshold which makes these 
states particularly interesting. Due to the suppression 
of the OZI-rule, these charmonium states usually have 
rather narrow widths. This makes their radiative tran- 
sitions and radiative decays having significant branching 
ratios and are experimentally accessible. It is also be- 
lieved to be the ideal hunting ground for exotic hadronic 
states like the glueballs whose existence is anticipated in 
QCD while its experimental signature remains obscure. 
Recently, the experimental interests have been revived 
with the upgrade for the BESIII experiment at BEPCII 
storage ring [l|, Q which collects charmonium samples 
that are orders of magnitude larger than ever. 

On the theoretical side, charmonium transitions have 
been studied using various methods. The physical pro- 
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cess involves both electromagnetic and strong interac- 
tions, the former being perturbative in nature while 
the latter being non-perturbative. Therefore, non- 
perturbative lattice calculations are preferred. Radiative 
transitions of charmonia have been studied comprehen- 
sively in quenched lattice QCD for the normal ground 
state charmonia Q and even for some excited and ex- 
otic ones [1]. However, an unquenched lattice study is 
still lacking. In this paper, we would like to pursue the 
feasibility of such a calculation using Nf = 2 dynamical 
twisted-mass fermion configurations generated by 

the European Twisted Mass Collaboration (ETMC). 

This paper is organized as follows: In Sec. [Hi we briefiy 
describe the lattice setup for the calculation of the hadron 
matrix element from three-point correlation functions in 
the theory. In Sec. IIIIl simulation details arc provided 
and the results are presented. This includes the charmo- 
nium spectrum and the dispersion relations, the single- 
quark form factors for rjc, XcO and the radiative transition 
matrix elements responsible for J/\l/ — ^ VcJ^ XcO ^ 
and he "Hcl- From these hadronic matrix elements that 
we obtained in our lattice calculation, we compute the 
transition decay width for these channels which are then 
compared with experimental values and the quenched re- 
sults. In Sec. II VI we will summarize our results and con- 
clude. 



II. THREE-POINT AND TWO-POINT 
CORRELATION FUNCTIONS 

The lattice setup in this calculation is analogous to the 
the vector form factor calculation of pions which has been 
studied extensively [T7l - l23| . Here we will briefiy review 
the general ideas involved. 

The transitions among charmonium states are trig- 
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gered by the electromagnetic interaction: Jl^lnt" 
J d^xAf^{x)jjl^'"^'\x) between the quark degrees of free- 
dom and the photon field. Here A^^{x) is the photon 

field and jjx '"^'\x) is the electromagnetic vector current 
of the quarks. Since the electromagnetic interaction is 
weak, one can treat it pcrturbatively. This leads to the 
computation of the hadronic matrix element of the cur- 
rent operator between the initial and the final ((/|) 

charmonium states: {f\jj:^'"^'\x)\i). We emphasize that 
although the electromagnetic interaction is perturbative, 



In this formula, interpolating operators which will cre- 
ate / annihilate the appropriate charmonium states are in- 
serted at time slices t = (the source operator) and t = t2 
(the sink operator), respectively. A local operator is used 
at the source while the sink operator with a definite three- 
momentum p2 is utilized. The current insertion at time 
slice t also carries a definite three-momentum q. Momen- 
tum conservation then implies that the initial state also 
has a definite momentum pi with q = P2— Pi- Physically 
speaking, the three-point function defined above repre- 



Thcrcforc. the desired hadronic matrix element 

{fiP2)\jj^'"^'\o)\'iiPi)) can be obtained once the 
energies Ei, E2 and the corresponding overlap matrix 

elements (il|02|/(p2)), (*(Pi)|Ci|f^) arc known, all of 
which can be obtained from corresponding two-point 



where Zi{p) = {^l\Oi\N{p)) is the corresponding overlap 
matrix element. 

With the relevant two-point and three-point functions, 

the hadronic matrix element (/(P2)|jl'^'™ ''(0)|j(pi)) 
could be extracted using two methods: The first is to fit 



the matrix element of the current between two hadronic 
states is in general non-perturbative. This is the quan- 
tity that we would like to compute using genuine non- 
perturbative methods like lattice QCD. 

Within the framework of lattice QCD, charmonium 
states arc realized by applying appropriate interpolating 
operators (Oi and O2 in the formula below) to the QCD 
vacuum \Q). Thus, the computation of the hadronic ma- 
trix clement {f\jif '"^'\x)\i) naturally leads to the follow- 
ing three-point function: 



(1) 

I 

sents a process in which an initial charmonium state with 
three-momentum pi created by 0| making a electro- 
magnetic transition to the final charmonium state with 
three-momentum p2 annihilated by O2 while the three- 
momentum difference q is carried away by the photon. 

Inserting a complete set of states between the electro- 
magnetic current operator and the charmonium opera- 
tors, one finds that, when t2 ^ t ^ 1, the states with 
the lowest energy dominate the three-point function: 



(2) 

I 

functions for the initial and final charmonium states. 

For this purpose, two-point correlation functions for 
the interpolating operators Oi for i = 1 , 2 are also com- 
puted in the simulation: 




I 

the two-point function Eq. ([3]) and three-point function 
Eq. ^ simultaneously. The second is to form an appro- 
priate ratio from the two-point and three-point functions 

and extract the matrix element (/(p2)|jl'^'™'''(0)|i(pi)) 
directly from the ratio. In this study, the second method 
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G^(<2,i;P2,Pi) - ^e-*P-''^e+''i--(f]|TO2(t2,X2)j('^-"-)(i,x)Ol(0,0)|r!) . 

X2 ,X 
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G,(t2,t;P2,Pl) \ (^^|02|/(P2))(»(Pl)|0ll»)(/(P2)U^""-HO)K(Pl)) • 

4£'l(.Plj£'2(,P2j 
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C,(t,p) ^^e-''P-(l]|O,(t,x)Ot(0,0)|r!) *^ ^M!e-^'(P)-? cosh 
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is utilized and the relevant ratio is defined as 



G,j(i2,i;P2,Pi 



C2(t2,P2) 
(/(P2)bf'"-^ 




t,Pl)C2ft,P2)g2ft2,P2) 

i, P2)Ci(i, pi)Ci(t2,Pi) 



4v/S2(P2)-Bi(pi) 



r 



(4) 



where the second line becomes valid when t2 3> t ^ 1, 
assuming only the corresponding ground states dominate. 
In this case, Rfj.it) becomes independent of t and fitting 
the ratio to a plateau behavior yields the desired hadronic 

matrix element (/(p2)|jM'' " ''(0)|i(pi)). 

Due to different implementations for fermions on the 
lattice, the electromagnetic current operator jji '"^'\x) 
might take different forms as compared with its con- 
tinuum counterpart. For Wilson-like fermions, like the 
twisted mass fermions that we use in this study, one could 
use either the local current or the conserved current. The 
local current is simpler in form but it is not conserved on 
the lattice. It thus requires an additional multiplicative 
renormalization given by the factor Zy, which of course 
can be determined non-perturbatively [3| . The conserved 
current is slightly more complicated but due to its con- 
servation, it does not need further renormalization, i.e. 
its multiplicative renormalization constant Zy = 1. In 
this work, we use the conserved current and the fact that 
Zy = 1 is also verified numerically in our simulation. 

In computing the three-point function defined in 
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Eq. ([2]), various quark contributions arise. Since the elec- 
tromagnetic current consists of contributions for all fla- 
vors of quarks, light flavors (i.e. u, d and s quarks) also 
contribute. Since our charmonium interpolating opera- 
tors are formed only from charm quarks, the contribu- 
tion from the light flavors can only occur through the so- 
called disconnected diagrams. The computation of these 
diagrams requires the light flavor quark propagators at 
basically all points on the lattice (the so-called all-to-all 
propagators). This is computationally extremely costly. 
Since the total electric charge of light quarks adds up to 
zero, one could argue that this contribution vanishes ex- 
actly in the flavor SU (3) limit. In this study, these contri- 
butions are neglected as is the case for previous quenched 
studies Q. Thus, we only need the charm quark contri- 
bution for the electromagnetic current which is propor- 
tional to the conserved current j^iix) on the lattice via: 

j'^'^""\x) = Qcjp,ix) with Qc being the electric charge 
of the charm quark. The conserved current jp,{x) for the 
twisted mass quark is given by 



j,i^) = cix)\-^U,{x)c{x + m) + c{x + ^)2li±lul{x)c{x) 



(5) 



When this current is inserted into the three-point func- 
tion, disconnected diagrams due to the charm quark in 
principle can still can arise. These are neglected in this 
study since charm quark is much heavier than the light 
quarks and they are also OZI-suppressed. Therefore, 
within the approximations described above, we only have 
to compute the connected diagrams from the charm cur- 
rent which can be treated using the sequential source 
method 



III. SIMULATION DETAILS 



The simulation setup for Nf 
fermions 



2 twisted mass 



Twisted mass fermions at the maximal twist are uti- 
lized in our study with two degenerate light flavors in the 
sea. The framework of maximally twisted mass fermions 
has been utilized in various studies of lattice QCD and 
are shown to be highly promising. It offers several advan- 
tages when tuned to maximal twist: (i) automatic 0{a) 
improvement [l^ is obtained when the bare untwisted 
quark mass is tuned to its critical value. Thus, only 
one parameter needs to be tuned; (ii) The determinant 
of the twisted mass Dirac operator is strictly positive. 
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protecting it against possible zero modes in the so-called 
exceptional configurations; (iii) It simplifies the operator 
mixing problem for renormalization. 

In this study, gauge field configurations using Nf = 
2 (u and d quark) twisted mass fermion are utilized. 
Other quark fiavors, namely the strange and charm 
quarks, are introduced as valence quarks. As discussed 
in refs. 0, [lol . fill [26j . we implement non-degenerate va- 
lence quarks in the twisted mass formulation by formally 
introducing a twisted doublet for each non-degenerate 
quark flavor. So, in the valence sector we introduce three 
twisted doublets, (u, d), (s, s') and (c, c') with masses /i;, 
/is and ^c, respectively. Within each doublet, the two 
valence quarks are regularized in the physical basis with 
Wilson parameters of opposite signs (r = — r' = 1). The 
fermion action for the valence sector reads: 

"5* = iXu,Xd) (Dw + merit + i/ii757"3) ( ) 

\XdJ 

+ {Xs,Xs'){Dw + merit +ifJ.s'-f5T3) { ] 

+ iXcXc) {Dw + merit + «Mc75T3) ( ] (6) 

KXc'J 

where D\y is the usual Wilson-Dirac operator, merit is 
the critical quark mass which is pre-determined from the 
simulation by ETMC. Numerically, it is more convenient 
to implement the operators in the twisted basis. For this 
purpose, one can perform a chiral twist and define the 
twisted doublets as follows: 

= <5-P(*^75r3/2)(^^) 

(J,) = exp (1^75X3/2) ^^''^^ 

(',) = exp{^u;J,r,/2)(^2) 

where lo — n /2 implements the full twist. 

In this work, all computations are done using Nf = 2 
twisted mass fermion configurations at the lattice spac- 
ing of a = 0.0666fm (/3 ~ 4.05). The size of the lattice is 
32'^ X 64 so that the spatial extent of the lattice is about 
2.13fm, which is a safe value for charmonium physics. In 
the temporal direction, anti-periodic boundary condition 
is applied for the quark field while periodic boundary 
condition is utilized in all spatial directions. The simu- 
lation parameters for our study arc summarized in Table 
III. 



TABLE I: Simulation parameters in this study. 



L -" *T 


/? 




a [fm] 




[MeV] 


-^conf 


32-" * 64 


4.05 


0.15701 


0.0666 


0.0080 


488 


201 



As for the charmonium states, we build the interpolat- 
ing field operators within the same Wilson parameters. 



In the physical basis, they read cTc and the correspond- 
ing form in twisted basis Xc^'Xc can also be obtained 
easily. These are tabulated in Table |ll] together with the 
possible J^*^ quantum numbers in the continuum and 
the names of the corresponding particle. ^ 

TABLE II: Local interpolating fields for charmonium states 
studied in this work in both physical and twisted basis, 
cTc = XcT'Xc- Also listed are the names of the corresponding 
particle and their J^'~" quantum numbers in the continuum. 









Xcl) 


Xcl 


he 


r 


li 


75 


1 


7i75 


a, J 


r' 


li 


1 


75 


7i75 


COi 




\'~ 


0-+ 


0++ 


1++ 


1+- 



Two-point functions are computed as usual for all char- 
monium states involved (those listed in Table |Il|) in our 
calculation. Fitting these two-point functions yields the 
energy for the corresponding charmonium states, both 
with and without three-momentum. As for the three- 
point functions, since only connected diagrams involving 
charm propagators are needed, sequential source method 
is utilized [2J]. The results for the two-point and three- 
point functions are then employed to construct the rel- 
evant ratio defined in Eq. (j4]). For definitcncss, we set 
t2 = 32 in our simulations which makes the three point 
function anti-symmetric (for jp=o) or symmetric (for ji 
with i = 1,2,3) [l3| about the time slice t2 = 32. In 
practice, we average the data from the two halves to im- 
prove statistics. All errors in this study are estimated 
using the conventional jack-knife method. 



B. Charmonium spectrum and dispersion relations 

Before computing the transition matrix element, the 
mass and the energy dispersion relations for the relevant 
charmonium states have to be verified. This is particu- 
larly important for our study due to the following rea- 
sons. Although charmonium spectrum has been studied 
extensively in quenched lattice QCD and the overall pic- 
ture agrees reasonably well with the experiment, some 
quantities like the mass splitting between rje and J/^' 
disagrees with the experimental value. It is widely be- 
lieved that this discrepancy mainly originates from the 
quenched approximation. It is therefore useful to check, 
using unquenched twisted mass configurations, whether 
this discrepancy can be resolved. Furthermore, although 
twisted mass configurations have been used successfully 
to study the light fiavors, using tliem on heavy charm 



^ In the first row of the table, we list the names of the charmonium 
states where XcO and Xcl are not to be confused with the charm 
quark field in the twisted basis. 
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quark needs some care. Being relatively heavy, the charm 
quark mass parameter //cO 0.2 in our study is not tiny. 
Of course, the good news from the maximally twisted 
mass lattice QCD is that it is 0{a) improved. Therefore, 
one would still hope to bring the lattice discretization 
errors under control. Another measure of the possible 
lattice artifacts is the charmonium mass, say the mass of 
the rjc meson m^^ in lattice unit. In our study, it turns 
out that rriri^a ^ 1. For charmonium states with non- 
zero three-momentum, this number becomes even larger. 
Therefore, one should carefully verify that these possible 
lattice artifacts for the charmonium states are not out 
of control. Only after these reassurances can one possi- 
bly proceed to calculate transitions among charmonium 
states reliably. As we will illustrate below, in our sim- 
ulation, most of these lattice artifacts are remedied by 
using the lattice dispersion relations for the charmonium 
states. 

Following Eq. ([3]), the energy E{p) for a particular 
charmonium state with three-momentum p can be ob- 
tained from the corresponding two-point function via 

The two point function is symmetric about t ~ T/2. 
In real simulation we average the data from two halves 
about t = T/2 to improve statistics. For each chan- 
nel, several three-momenta (including the zero three- 
momentum) have been computed. Different momentum 
modes that are related by lattice symmetries are aver- 
aged over. 

The effective mass plateaus at zero three-momentum 
for the charmonium states studied in this work are il- 
lustrated in Fig. [TJ From lower to higher values, the 
plateau corresponds to the charmonium state of rjc, J/^ , 
XcOj he and Xci, respectively. It is seen that the effec- 
tive mass values for rjc and J/^" have shown very clear 
and well-established plateau behavior, resulting in rather 
small statistical errors. We use the mass of J/^* from our 
simulation to set the bare charm quark mass parameter 
fic- After some tuning, we find a/ic — 0.203 roughly cor- 
responds to the value that is consistent with the value 
quoted in the Particle Data Group (PDG). We then fix 
12c at this particular value for all our subsequent calcula- 
tions. Since twisted mass lattice QCD is 0{a) improved, 
the anticipated cutoff effects induced by the charm quark 
mass is roughly ©(a^/x^), which is at a few percent level. 
This of course still needs further verification from mea- 
sured physical quantities. The effective mass plateaus for 
other charmonium states: XcOi Xci £md he are relatively 
noisy with larger statistical errors. The fitted effective 
mass values are collected in Table IIIII which can be com- 
pared with the corresponding values from PDG. 

It is gratifying to see that our lattice result suggests 
rrij/^ — nir]^ = 104MeV, a little smaller than the PDG 
value of about 117 MeV. This is already a great improve- 
ment over the quenched studies where the lattice results 
are typically away by dozens of McV. This remaining dis- 



TABLE III: Charmonium effective mass [Unit:MeV] 









XcO 


Xci 


he 


Mass (this work) 


2997.4 


3101.0 


3326.9 


3613.2 


3466.8 


Error 


0.5 


0.7 


4.5 


18.4 


23.0 


PDG 


2980.3 


3096.9 


3414.7 


3510.7 


3525.9 



crepancy might come from lattice artifacts (since we are 
simulating at a fixed lattice spacing without taking the 
continuum extrapolation) and/or from the fact that we 
have neglected annihilation diagrams for the charm quark 
in the two-point function, as estimated in Ref. pTj . 

To get a feeling about the size of the lattice artifacts 
for the charmonium states with non-vanishing three- 
momenta, we investigate the dispersion relations for 77c, 
J/^* and XcO states. The energy i?(p) are obtained from 
the corresponding effective mass plateaus of the two- 
point functions with prescribed three-momentum. As 
said in the beginning of this subsection, since the char- 
monium states are relatively heavy in lattice unit, the 
continuum dispersion relation E"^ = m? -\- c^p^ may not 
be a good description, where c is the speed of light which 
should be close to unity if lattice artifacts are small. In- 
deed, our data suggest that the naive continuum disper- 
sion relation is violated with the fitted value of sub- 
stantially away from unity by as much as 12% even for 
•qc and J /'^ states. However, we find that, if we utilize 
the standard lattice dispersion relation 

4sinh^(^)^4sinh^(f)+Zx45:sin^(|) , 

(9) 

which recovers the naive dispersion relation in the con- 
tinuum limit, we could describe our data extremely well 
with the fitted values of Z for rjc and J /'^ rather close to 
unity. 

The dispersion relations for t^c, J/^* and XcO are illus- 
trated in Fig. [21 Fig. [3] and Fig. HI respectively. We find 
that Z^^ = 1.063(8), Zj/vt = 1.056(7), Z^^^ = 1.13(24), 
all of which are close to the anticipated value Z = \. The 
difference seems to be at the order of 0{{jico)^) ^ 4% as 
the naive estimate suggests. In evaluating the two-point 
function, thanks to the averaging of various p related by 
lattice symmetries, we get very good dispersion relation 
even at p^ as large as 10p^j„, where Pmm = (100) (in 
unit of {2tt) / L) is the minimal lattice momentum. It is 
also seen that even at the largest three-momentum, the 
lattice dispersion relation still offers a very good descrip- 
tion of the data. This gives us confidence that, at this 
particular lattice spacing that we are simulating, most of 
the lattice artifacts for the charmonium states are taken 
care of by using the lattice dispersion relation 

C. Form factors for rjc and XcO 

In the continuum, the hadronic matrix element 
(^c(P2)|i/j(0)|?7c(Pi)) may be parameterized by only one 
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1.3 
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CO 



1.1 
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Tic ^ — I 1 

XcO 
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FIG. 1: Charmonium effective mass plateaus. From lower to higher values, the plateau corresponds to the charmonium state 
77c, J/5', XcO, he and Xci, respectively 



Ti;, Dispersion Relation 




0.15 0.2 0.25 

4*Ii sin^(Pi /2) 



FIG. 2: The t^c dispersion relation obtained from our calculation. Following the lattice dispersion relation ((9|, the quantity 
4sinh^ (""2^) (vertical axis) is plotted versus different values oi A^^sin^ (^) (horizontal axis). The data points with errors 
are simulation results while the straight line is a linear fit according to Eq. ((5)) with the fitted value of Z indicated in the upper 
right corner of the plot. 
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J/y Dispersion Relation 




4*Ii sin^(Pi 12) 



FIG. 3: The same as Fig.[2]but for J/*. 



form factor f{Q^) as: Q 

{vc{P2)\j^.mvc{Pi)) = /(Q')bi +P2)m , (10) 

where = — (p2 ^Pi)^ is the square of four momentum 
transfer. This quantity is also called the single-quark 
elastic form factor in Ref. [28|. It is not a directly mea- 



surable quantity experimentally. But it is a quantity that 
can be computed in lattice simulations which can then be 
utilized to compare with similar results from models (see 
Ref. d^D- Note that this form automatically ensures 
the current conservation (??c(P2)|c^^.7p|'7c(Pi)) = since 
Q- iPi+ P2) = ij>2 - Pi) ■ iP2 + Pi) =0. On a finite lat- 
tice, the partial derivatives are replaced by correspond- 
ing finite differences on the lattice. For the temporal 
components of the four-momenta, this amounts to re- 
placing the continuum energy by its lattice counter part: 
{Pi)o 2 sinh(i5j;/2). Note that this modification applies 
to the energy factors inside the square root in the second 
line of Eq. @. Of course, in principle the spatial com- 
ponents should also be modified according to the lattice 
dispersion relation ^ . But since our three- momenta are 
relatively small in lattice unit, this replacement does not 
make a significant change. For the temporal components, 
however, since aE{p) ^ 1 for all charmonium states be- 
ing studied, this modification is crucial. For example, 
according to Eq. (fTO)) . the form factor f{Q^) is a scalar 
function which is the same for all indices fi ~ 0, 1,2,3. 
Only after using the modifications suggested by the lat- 
tice dispersion relation can we obtain consistent results 
for f{Q^) at different values /i. 

To obtain the desired hadronic matrix element 
('yc(P2)|jp(0)|77c(pi)), we form the ratio defined in 
Eq. This is done for the zero three-momentum 

case P2 = (0,0,0) and for various non-vanishing three- 
momenta. In Fig. [5] and Fig. [BJ we display the typical 
behaviors for Ro{t) for p2 = (0,0,0) and p2 = (0,0,1), 
respectively. It is seen that clear plateau behaviors have 



been established from which the form factor f{Q^) can 
be extracted. We have checked that taking the tempo- 
ral and spatial components of the current yields consis- 
tent results for f{Q^) although the results obtained from 
11 = (i.e. Ro{t)) gives smaller statistical errors, which 
we take as the final result for the form factor at that 
particular Q^. 

The fitted values of f{Q^) obtained from the ratio are 
shown in Fig. [7] versus different values of where two 
different type of symbols stands for p2 = (0,0,0) and 
P2 = (0,0,1), respectively. It is seen that the data ob- 
tained in the two cases tend to lie on a universal curve. 
Following Ref. Q , we fit the data for the form factor with 
the following function: 



/(g2) = cxp 



16/32 



(1 + aQ' 



The fitted parameters turn out to be: 

a = -0.096(6) GeV"^ ^ = 567(2) MeV 



(11) 



(12) 



This value of /3 is larger than the corresponding value 
480(3)MeV obtained in the quenched approximation in 
Ref. [3| , making the corresponding form factor obtained 
from our unquenched calculation "harder" (i.e. decays 
slower with increasing Q"). The comparison of this form 
factor obtained from various phenomenological models 
with the corresponding quenched result has been ad- 
dressed in Ref. 0, [l^ . It was noted that using the simple 
harmonic oscillator (SHO) wavefunctions yields a harder 
form factor when compared with the quenched lattice 
result. In the quark model of ISGW [23, however, an 
extra factor k, ~ 0.7 was introduced such that the form 
factor takes the form f{Q^) ~ exp (— (5'^/(16/3^k^)) near 
= which agrees with the quenched lattice result well, 
given a phenomenological value of /3 ~ 710MeV. Since 
our unquenched lattice result suggests a harder behavior 
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FIG. 4: The same as Fig.[2]but for XcO- 



P2=000 
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FIG. 5: The ratio Ro{t) defined in Eq. (g]) for 77^ with p2 = (0,0,0). 
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FIG. 6: The ratio Ro{t) defined in Eq. (g]) for 77c witfi p2 = (0,0, 1). 



for the form factor than the quenehed case, we find that, 
for the same value of (3 taken in the model, a factor of 
K ~ 0.8 will make the model predictions in good agree- 
ment with our unquenched lattice results. One can define 
a squared mean charge radius (r^) with (r^) given by: 



-6 



I6/32 



(13) 



Our unquenched lattice result then yields (r^) ^ 
0.213(1) fm which is smaller than the corresponding 
quenched value of G.255(2)fm. 

Our unquenched result yields a harder behavior for 
the form factors which can be understood qualitatively. 
Physical scales on the lattice are usually set by some 
long-distance physical quantities, like the static quark 
anti-quark potential in the quenched or the pion decay 
constant in the Nf = 2 twisted mass lattice QCD. How- 
ever, it is known that quenched lattice QCD did not 
reproduce the true QCD /3-function due to the lack of 
the quark loops. In particular, when running from the 
lower energy scale up to the scale of charmonium physics, 
quenched lattice QCD gives a weakened strong coupling 
constant than unquenched lattice QCD. This is believed 
to be the major reason for the discrepancy between the 
mass splitting of J/^* and rjc in quenched lattice QCD 
with the true experimental result. Therefore, unquench- 
ing the quarks will basically make the effective coupling 
constant stronger at charmonium scale when compared 
with the quenched case. This in turn gives a smaller 
charge radius for the unquenched case, in agreement with 
what we find in our calculation. 

The hadronic matrix element (Xco(P2)| jM(0)IXco(pi)) 
for Xco has the same form of decomposition as that for rjc- 
The corresponding form factor is defined as in Eq. ([TU]). 
In exactly the same manner, we can obtain the form fac- 
tor f{Q^) for XcO except that we have only computed the 



case p2 = (0, 0, 0). This is illustrated in Fig.[8l The data 
is fitted with the function: 



.f(Q') = .f(0)exp 
The fit parameters are: 



I6/32 



/(O) 1.0002(5), /3 = 510(16) MeV 



(14) 



(15) 



This value of /? is also larger than the quenched value 
of 393(12)MeV from Ref. [J] (i.e. the unquenched form 
factor is also harder than the quench one). Note that 
our fitted value of /? for x^o is smaller than that for 77^ 
resulting in a larger charge radius for XcO when compared 
with that for 77c. This is consistent with the quark model 
picture since XcO in this model is a i = 1 state and the 
charge radius is naturally larger due to the presence of 
the centrifugal potential. 

Finally, we remark that, in cases of both rj^ and XcOj 
although some data might have large errors at non- 
vanishing Q^, the form factor as zero momentum transfer 
is always consistent with unity which is in fact a man- 
ifestation of the current conservation: Zy = 1- There- 
fore, the conserved current requires no extra multiplica- 
tive renormalization as it should. This is verified numer- 
ically by our simulation data. 



D. J/'I' — s- r]c"f transition 

The matrix element (7/c(P2)|j^(0)| [J/^']^ (pi)) is re- 
sponsible for the calculation of J/^* — ijcj transition 
rate. Here we use the index r to designate the polar- 
ization of the initial J/'I' state whose polarization vector 
is denoted by e^(pi,r). In the continuum, this matrix 
element can be decomposed as Q 
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FIG. 7: The form factor f{Q^) for rjc obtained from p2 = (0,0,0) (red data points) and p2 — (0,0, 1) (blue data points). The 
curve is a fit for all the data using the functional form of Eq. (|11[). 
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FIG. 8: The same as Fig. [7| except that the charmonium state is XcO- 



{r,M\.r m [J Mr ipi))^ 



(16) 



Thus the matrix element is characterized by one form 
factor V{Q^). By forming the appropriate ratio, relevant 
lattice results V{Q^) are extracted from the plateaus of 
the ratios. The relation of V{Q^) with its continuum 
counterpart V{Q'^) is V{Q^) = 2 x |e x V{Q'^), where 
the factor 2 comes from the quark and the anti-quark 
while the factor (2e/3) is due to the charge of the charm 
quark. The results for the transition form factor V{Q^) 



thus obtained are illustrated in Fig. [3] Following Ref. 
the data is fitted with the function: 



16/32 



(17) 



The resulting fitted parameters we find are as follows: 
y(0) = -2.01(2), /3 = 580(19) MeV (18) 
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This is to be compared with similar results from pre- 
vious quenched lattice study: V^(0) = —1.85(4) and 
/3 = 540(10) MeV in Ref. [|. 

With the values of the transition form factor on the 
lattice, the J/^E" — >■ rjc^ decay width can be obtained: 



r (J/* ^ 77c7) = a 



64 



27 (m^^ + m^)^ 



1^^(0)1 



(19) 



where q is frame dependent. If we choose the frame 
in which the initial J/^* is at rest, we have: |qp ~ 



( 



mz. 



/(4m|,). Substitute this into Eq. ([T9]) wc 



then get the J/'^ rjd decay width: 



2.84(6) KeV, T„ 



1.99(6) KeV (20) 



where ,^ denotes the result with physical mass values 



(e.g. values from PDG) are substituted into Eq. (fT9l) . 
while Fmiat stands for using the mass values computed 
from the lattice directly. This difference arises since our 
lattice results for the masses for J/^* and rjc in subsec- 
tion IIIIBI (see Table IIIip do not coincide with their ex- 
perimental values exactly. Although our lattice value for 
the mass of the J/^P is quite close and that for r/c is also 
closer to the experimental value than the corresponding 



quenched value, the decay width turns out to be propor- 
tional to (rri j/^ — rrijj^)^ which magnifies the difference. 
Note that, for quenched lattice calculations, using dif- 
ferent charmonium mass values makes a even bigger dif- 
ference, as noted in Ref. The corresponding results 
are: E™^,, = 2.57(11) KeV, E™,^. = 1.61(7) KeV. This 
is due to the fact that quenched lattice calculations yield 
a much smaller value for mj/q, — rrir^^ when compared 
with the true experimental value. In our unquenched 
study, however, we see that this difference is somewhat 
milder compared with the previous quenched situation. 
Both our lattice result and the previous quenched re- 
sult for this quantity are to be compared with the value 
TpDG = 1.58(38) KeV quoted by the PDG. Note that the 
PDG value is an average of CLEO result and the Crys- 
tal Ball result, the former being 1.92(30) KeV which is 
closer to our lattice result while the latter from Crystal 
Ball being 1.18(33) KeV, smaller than lattice results. 



E. XcO — !• J/*I'7 transition 

In the continuum, this transition matrix element has 
the following decomposition: 0] 



r,2 



+ '^"'Z-' mve{pv,r) ■ ps[pv ■ Ps{pv +PsT - m^sPv - ^vPs\ 



(21) 



with ^l{Q'^) = {pv-ps) — myrrig. Therefore, the 
hadronic matrix element is characterized by two form fac- 
tors Ei{Q'^) and Ci{Q^). At the physical photon point 
with = 0, only the former contributes. 

The form factor Ei{Q^) can be obtained by following a 
similar process as the other form factors. We can always 
choose some combinations of pv , Ps such that 

{S{vsWmV{vv,r))<xE,{Q^). 

The final lattice results for Ei{Q^) are shown in Eig. [TOl 
Wc then use the following form 



Ei{Q') = Ei{Q)\l + %\e^V 
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(22) 



to fit the data [3|. The fitted parameters we obtain are: 

aEi{Q) = -0.1699(51) , (23) 
p = 871(85) McV, P = 451(62) MeV . 

The fitted value of i?i(0) at = Q is also indi- 
cated in Eig. [TU] together with its error. These results 



are to be compared with similar results from previ- 
ous quenched lattice study 3: atEi{Q) ~ —0.137(12), 
/? = 542(35)MeV and p = 1.08(13) GeV. 

At the physical photon point ~ 0, the decay width 
for this radiative transition is given by: 



E(xcO ^ J/*7) - a^^|^i(0)|2 
9 mi 



(24) 



where Ei is related to Ei by: 



^i(Q') = 2x -ex^i(Q2) 

Substituting our lattice result for i?i(0), we then can get 
the decay width in physical unit: 



= 85(7) KeV, E,„,^, 



65(4) KeV (25) 



which is to be compared with the quenched lattice result 
of E„pj,^ = 232(41) KeV and E„,^^ = 288(60) KeV. It 
is seen that our unquenched result for this decay width 
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O -1 




Q''(GeVT 

FIG. 9: The lattice results for the J/^ — >■ 77c7 transition form factor V{Q^). The curve is a fit using the function in Eq. (|17p . 
The fitted value of V{Q^ = 0) is also shown at = together with its corresponding error. 



is substantially smaller than their quenched values. The 
value for this quantity quoted by the PDG is given by: 
TpDG = 119(11) KeV which lies in between the quenched 
and unquenched results. 



he 



J7c7 transition 



The form factor decomposition for this process is iden- 
tical to Eq. ([2T|) . However, the signal for the state he is 
much noisier. It turns out that we could only get rea- 
sonable signal with ~ 000. In order to get various 
values of Q^, we vary the values of q and p,,^ simultane- 
ously such that he is always at rest. 

The form factor we obtain is illustrated in Fig. [TlJ We 
fit the data with a functional form: 



16;32 



(26) 



and the fitted parameters come out to be: 

aEi{0) = -0.39(1), 13 = 440(23) MeV (27) 

The fitted value of Ei (0) is also shown in Fig. [TT]at = 
together with the corresponding error. These numbers 
are to be compared with the corresponding quenched re- 
sults: atEi{0) = -0.306(14) and ^ = 689(133) MeV 
The physical decay width for the transition is given by: 



r {he ricj) ^ a 



16jq[ 
27 < 



1^1 (0)1 



(28) 



With our lattice result for £"1(0) substituted in, we find: 



''■phy 



234(12) KeV, F^ 



210(13) KeV . (29) 



The corresponding quenched lattice values are: Fm^j^^ = 
601(55) KeV and F^^^^ 663(132) KeV, both of which 
arc about a factor of 3 larger than our unquenched re- 
sult, though the errors are somewhat large. The lattice 
results for this decay can now be compared with the re- 
cent measurement at BES-III. 0^ The total width for he 
and the corresponding branching ratio for the radiative 
transition are found to be 

F*°'p = 730 ± 450 ± 280 KeV , 
B{he^rie-/) = (54.3 ± 6.7 ± 5.2)% , 

where the first error is statistical and the second is sys- 
tematic. If we multiply the central values for the above 
two quantities and add the errors in quadrature we find 
the decay width r{he Tjel) = 396(294)keV, which 
could be compared with our lattice result. The agreement 
within a large error is seen although improvements from 
both experiment and lattice calculations are required to 
cut down the large uncertainties for this quantity. 



IV. SUMMARY AND CONCLUSIONS 

In this exploratory study, we calculate the form factors 
for some of the ground state charmonia and their radia- 
tive transitions using unquenched Nf = 2 twisted mass 
fermions. The s and c quarks are quenched which arc in- 
corporated via a twisted doublet for each non-degenerate 
quark flavor in the valence sector. Our study focuses 
on the form factors for rje, XcO and the J/^P — >■ r]e"f, 
Xco ^ J 1^1, he — !■ 77c7 radiative transitions. The mass 
spectrum and dispersion relations for these charmonium 
states are first examined. Good agreement of the com- 
puted spectrum with the experiment is found. It is also 
verified that, by using lattice dispersion relations instead 
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TABLE IV: Summary of the results obtained in this work. Also listed are the corresponding results from quenched lattice 
QCD 0]. Experimental values or values from PDG are also listed whenever available. 



Fitted parameter /3[MeV] for form factors 





r?c 




This work 


567(2) 


510(16) 


Ref. [3| 


480(3) 


393(12) 



/Fmiat [keV] for transitions 





J/* Tjc-y 


Xco J/*! 


he lid 


PDG 
This work 
Ref. [3] 


1.58(38) 

2.84(6)/1.99(6) 

2.57(11)/1.61(7) 


119(11) 

85(7)/65(4) 

232(41)/288(60) 


396(294) 

234(12)/210(13) 

601(55)/663(132) 



of the naive continuum ones, the lattice artifacts for these 
charmonium states are well under control. By comput- 
ing various appropriate ratios of the three-point functions 
to the two-point functions, hadronic matrix elements for 
these transitions are obtained at various of . Using 
the parameterized form in terms of relevant form factors, 
we obtain the lattice results for the relevant form factors 
and the radiative decay widths for these channels. Our 
results are summarized in Table HV] which are to be com- 
pared with those obtained in previous quenched lattice 
studies and experimental values. 

Although some quantities from our unquenched study 
turn out to be comparable with the quenched results, 
quite a number of the results still differ substantially, as 
is seen from Table IIVI For example, for the form fac- 
tors of 7?c and XcO a harder behavior (a larger value of 
j3 hence a smaller charge radius) than the quenched re- 
sult is found. As for the decay width for J/'I' — > r^cT, a 
value larger than the quenched result is obtained. Due 
to the improvement of the mass splitting for the char- 
monium in unquenched study, the discrepancy between 
the results using the physical mass and the lattice com- 
puted mass is somewhat narrowed, although the value is 
still larger than the value quoted by PDG. For the decay 
width of XcO ^ J/'^l, a value smaller than the PDG 
value (and also the quenched result) is obtained. As for 
the he — >■ 77c7 transition, the signal is noisy and our un- 
quenched result is much smaller than the quenched value 
with a large statistical uncertainty. It is still compati- 
ble with the recently measured value at BES-III which 
also has a large error. To get better signals for this chan- 
nel, variational methods or smearing techniques might be 
necessary which will be investigated in the future. 

In this preliminary study, we simulate at only one lat- 
tice spacing and sea quark mass, and no chiral nor con- 



tinuum extrapolation is made. The physics involved in 
this study mainly concerns the heavy flavor part of the 
theory which should not be sensitive to the pion mass. As 
for the lattice artifacts, we argued that, thanks to the au- 
tomatic 0{a) improvement, the lattice artifacts is under 
control. Indeed, by using the lattice dispersion relations, 
we verified that all charmonium states that we studied ex- 
hibit controlled lattice errors in their dispersion relations 
of about a few percent, which is roughly at the order of 
(/ico)^ for our simulation. With the experience gained in 
this study, it would be better and also possible to study 
charmonium radiative transitions in a more systematic 
manner (more lattice spacings, more pion mass values 
etc.) using unquenched lattice QCD. It is also tempting 
to perform similar studies with the Nf = 2-|-l-|-l dynam- 
ical twisted mass fermion. Given the promising experi- 
mental status of BES-III at BEPC-II, the unquenched 
lattice studies on charmonium transitions will certainly 
be an interesting project to pursue in the future. 
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